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ABSTRACT 

Red-sequence galaxies record the history of terminated star-formation in the Universe and 
can thus provide important clues to the mechanisms responsible for this termination. We con- 
struct composite samples of published cluster and field galaxy photometry in order to study 
the build-up of galaxies on the red-sequence, as parameterised by the dwarf-to-giant ratio 
(DGR). We find that the DGR in clusters is higher than that of the field at all redshifts, im- 
plying that the faint end of the red-sequence was established first in clusters. We find that the 
DGR evolves with redshift for both samples, consistent with the "down-sizing" picture of star 
formation. We examine the predictions of semi-analytic models for the DGR and find that 
neither the magnitude of its environmental dependence nor its evolution is correctly predicted 
in the models. Red-sequence DGRs are consistently too high in the models, the most likely 
explanation being that the strangulation mechanism used to remove hot gas from satellite 
galaxies is too efficient. Finally we present a simple toy model including a threshold mass, 
below which galaxies are not strangled, and show that this can predict the observed evolution 
of the field DGR. 
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• 1 INTRODUCTION 

' Red-sequence galaxies are an important population for understand- 
ing galaxy formation and evolution in general. They represent sys- 
tems with very little or no on-going star-formation and thus are a 
unique tracer of the past activity of galaxies. De Lucia et alj | |2004) 

, found a deficit of faint red-sequence galaxies in clusters at z~0.8 
relative to local clusters. This would suggest that star formation 
ended earlier for the most luminous/massive galaxies at high red- 
shi ft and would supp ort the "down-sizing" picture first proposed 
bv lCowie et alj ( Il996^ . in which the termination of star-formation 
progresses from the most massive to the least massive galaxies as 
the universe ages. Similar results for clusters spanning a range of 
redshifts have been found by Tanaka et al. (2005), de Lucia et al. 
l2007l) . IStott el^alj fcoOTh . loilbank et alj HOOT,) and others. Simi- 
larly, surveys of field galaxies have attempted to characterise the 
evolution of the red-sequence luminosity function. The relative 
contributions of passive evolution/termination of star-formation 
and number density evolution/dry mergin g are still open questions 
(e.g.. lBell etaLll2004lcrmatti et al.ll2006l) . 

Red-sequence galaxies can thus place constraints on galaxy 
formation prescriptions in theoretical model s. Recent observa- 
tions in the local univ erse using SDSS data l iBaldrv et alj|2006l : 
IWeinmann et al J [20061) have found that the fractions of satellite 
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galaxies which are red are too high in current semi-analytic models. 
With the recent abundance of surveys targeting red-sequence galax- 
ies both in the cluster and the field, the time is now ripe to draw 
together these works to build a picture of the build-up of galaxies 
on the red-sequence as a function of environment and to confront 
these with predictions of galaxy formation models. 



2 OBSERVATIONAL DATA 

We aim to construct the most uniform sample possible using 
colour-selected red-sequence galaxy data from the literature. To 
this end, we transform all our data onto the same system. A simpler 
quantity than the luminosity function (LF), which has been em- 
ployed in the past, is the red-sequence dwarf-to-giant ratio, here- 
after DGR. This is essentially just a LF reduced to two bins and 
avoids the complications of having to fit an analytic function (usu- 
ally with d egeneracies between t he fitted parameters) to the data. 
Following jPe Lucia et alj j2004h . we define luminous, or giant, 
galaxies as those brighter than Mv ~ —20.0 and faint, or dwarfs, 
as those with —20 < My ^ —18.2 (k- and evolution-corrected). 
Wh ere data are p resented on a different photometric system, we 
use iBruzual & C hariot (2003) stellar population models to calcu- 
late the transformation using a Zf = 3, Solar metallicity, single 
stellar population with a Chabrier IMF. Where DGRs have been 
given in the original paper using these limits (as is the case for 
much of the cluster data) we use these numbers and their associated 
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errors directly. For other works, we take the M* and a values from 
the Schechter function fit, transform to My and integrate the trans- 
formed LF to calculate DGR. We estimate conservative errors using 
a Monte Carlo method taking the errors in the fitted parameters and 
allowing for their covariance by employing the error ellipse given 
in the original reference. In addition we propagate a 0.05 magni- 
tude uncertainty (discussed below) in the magnitude limits to allow 
for any residual mismatch in the converted photometric system. 

For the field samples, we carefully examine all of the 
Schechter function fits in each redshift bin and only use measure- 
ments wher e the data ful l y cov er the dwarf and giant range. For 
example, the lBrown eTallllOOTl NOAODWFS) data do a good job 
of constraining the shape of the LF over the magnitude range we 
require for but at higher redshift are in sufficiently deep to 

constrain the number of dwarfs. Conversely, the lZucca et al.l ( l2006l 
VVDS) data are sufficiently deep to measure the faint end of the 
LF at higher redshift, but at z^O.6 are subject to considerable un- 
certainty around A/*, presumably due to the relatively limited area 
of the survey . The r emainder of our field samp le uses p oints from 
iBaldry et all ( |2004 SPSS). ISell et alj 1 I2OO3I , SDSS), ISell et al.l 
(I2OO4I COMB O-17).lDriyer et al.l (l200d Millennium Galaxy Cata- 
logue) and Sca rlata et al.l ( l200l COSMOS). Fig. 1 shows DGR as a 
function of redshift for these surveys. Where the redshift ranges of 
different surveys overlap, we find good agreement in our measured 
DGRs (e.g., VVDS and COSMOS a.t z = 0.7). 

Clusters show c onsiderable scatter in their individual DGRs 
(e.g.. Ide Lucia et al] i2007) and so we only focus on works mea- 
surin g DGR from ensemb l e averages of clusters in each redshift 
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Two results are immediately apparent from the observational 
data in Fig. 1: the DGR of field red-sequence galaxies is always 
lower than that of cluster galaxies at the same redshift, and the 
DGR of both samples evolves with redshift. We parameterise the 
evolution of the DGR with a fit of the form (1 + z)^ finding 
/3 = -1.8 ± 0.5 and (3 = -3.4 ± 1.1 for the cluster and field 
samples respectively. 

2.1 Sources of systematic error 

The first systematic to be addressed is that of transforming pub- 
lished photometry on different systems to a uniform one. Most of 
the published cluste r studies give rest-frame V magnitudes (e.g., 
de Lucia et al. 2007), and most field studies give rest-frame B (e.g.. 
Bell et al.. 2004 ). This is not exclusively the case, however, and we 
find good agreement where different studies measured in differ- 
ent pass bands in the same environments at similar redshifts over- 
lap (e.g.. [Hansen et a l. 2007; Stott et al. 2007). Thus we believe the 
transformation is not a large source of error. We have propagated 
a generous 0.05 magnitude error through our analysis to allow for 
this. This also allows for the dominant transformation uncertainty 
which is that of the systematic calibration in the individual sur- 
veys (e.g., Hansen et al. 2007 quote their calibration un certainty as 
0.03 m ag). In fact, comparing Hansen et alj |2007|) and lStott et al.l 
1 I2OO7I) , the measured DGRs agree within the errors even if this ad- 
ditional uncertainty is neglected. So, we believe that our error bars 
are ample. 

The important other systematic for field studies is cosmic vari- 
ance. The authors of the various works included here discuss the 
size of this effect for their surveys, and this is not likely to be a 
major contribution to the error bars in Fig. I. Cluster studies re- 
quire more care in the sample selection to minimise systematic 



errors. As found by e.g., iGilbank et al] tOOH) and [Hansen et al.l 
1 I2OO7I) , the faint end slope of the cluster RSLF depends on cluster 
richness/mass (at least for z<0.5), in the sense that low mass clus- 
ters/rich groups are expected to have systematically lower DGRs. 
In this compilation, we have attempted to utilise the richest/most 
massive systems at each redshift for the cluster sample, to maximise 
the difference between the clusters and the field. In so much as, on 
average, the most massive systems at one redshift evolve into the 
most massive systems at lower redshift, such a selection is appro- 
priate for evolutionary studies. If we consider the three data points 
nearest to z~0.5, these were drawn from very differently selected 
samples: a statistical sample of 54 red-sequence selected clusters, 
having richnesses co rresponding to a mass limit of ;> 6 X lO^M© 
jGilbank et alj2007l) ; a sample of 1 1 clusters selected via their un- 
resolved clusterlightfroni_sh optical imaging, having various 
masses jde Lucia et al. 2007h: an d a sample of ten of the most X- 
ray luminous clusters jStott et aT 2007). The good agreement of the 
points from these ensemble clusters suggests that the DGR is rela- 
tively insensitive to mass or selection method over this mass range. 
Furthermore, it suggests that the evolution of the cluster DGR is 
likely not due to a systematically changing mass limit with redshift. 

Additionally, the virial radi us of a cluster depend s on its mass 
roughly M^^''^ (e.g.. ICarlberg et al.iri997h . so using a 

fixed physical radius within which to measure a property of a clus- 
ter probes a different fraction of the cluster's extent depending on 
the mass of the cluste r. The effect of this on DGR is small as can 
be seen from fig. 5 of Ide Lucia et al.l ^007). By limiting the mass 
range of the clusters to the most massive systems, we have min- 
imised the possible impact of such an effect, since the difference 
between a fixed physical radius and a fixed fraction of the virial 
radius only depends on the mass to the 1 /3 power. 

2.2 Conflicting results 

lAndreo j JioOTh recently argued, based on a small number of indi- 
vidual clusters, that the cluster DGR does not appear to evolve with 
redshift. We indicate his data points on Fig. 1, omitting the three 
z> 1 clusters for which the data do not bracket the 4000A break. 
All of the individual clusters are in reasonable agreement with the 
ensemble averages, and with our best fit line, except for one at 
z = 0.89. Furthermore, his sample only contains two z < 0.5 clus- 
ters, one of w hich is Coma, f o r whi ch the DGR is in good agree- 
ment with the Ide Lucia et al. I ( l2007h value. Given the significant 
cluster-to-cluster scatter, the ensemble samples compiled here pro- 
vide a better reflection of the global average in clusters, and allow 

us t o measure sign i ficant evolut ion with redshif t. 

iTanaka et al] ( |2005[) and Ide Lucia et al] ilOOH) both con- 
structed composite s amples of local S DSS clust ers from the 
C4 clu ster catalogue jMiller et al]|2005l) . However, iTanaka et al.l 
l l2005h found a low DGR, comparable with their two higher red- 
shift cluster z=0.55 and z=0.83) and therefore concluded that the 
cluster DGR might not evolve with redshift. Both of their higher-z 
individual clusters lie on the best-fit r elation of Fig . 1, so the dis- 
crepancy arises from the Z'^O point, ide Lucia et ahi ( l2007t) found 
a higher DGR (for th eir a > 600 km s~ ^ sample) more consis- 
tent with the results of lHansenet"ai] ( l2007l) for the richest clusters; 
IStott et al.l koot for the most X-ray luminous clusters; and also 
for the regions of highest density in the SDSS (Baldry et al. 200^ 
(open sq uare). We sug gest that the reason for the low DGR mea- 



sured by ITanaka et alj l l2005l) may be that their sample was dom- 
inated by lower mass systems. We note that their measurement 
is just consistent with the lowest DGR point at Z'^O. 1 in Fig. 1 
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Figure 1. The evolution of the Dwarf-to-Giant ratio, DGR, for vaiious samples taken from the literature and converted to a uniform system, as described 
in the text. Fille d circles and open diamonds denote cluster (Barkhouse et al. 2007, de Lucia et al. 2007, Gilbank et al. 2007, Hansen et al. 2007, Stott et al] 
l2007t) and field iBell et alj|20M.lBaldrv et alj|2004 iBell et al.l i2004. Driver et al. 2006. ,Zucca et al.ll200a IBrown et aUl2007i. iScarlata et alj|2007l) samples 
respectively from observational data, and lines show (1 + z)^ fits. Horizontal error bars indicate the redshift range covered by the data. Naked error bars are 
measurements for individual clusters from Andreon (2007) and jTanaka et al J ^200^. labelled 'A' an d 'T' respectively. Only ensemble clusters are included 
in the fit. The filled square shows t he highest density r egions in the SDSS from lBaldrv et alj j2006l) which agrees well with the cluster data. Asterisks and 
triangles show predictions from the lBower et al.l i200d) semi-analytic model, connected with lines to guide the eye. The dot-dashed lines show the predicted 
evolution of the field DGR if isolated galaxies have a constant DGR=0 and the observed trend is due to the increasing abundance of groups (having the same 
DGR as clusters) above the labelled threshold mass, Mth with cosmic time, as described in §4. 



which comes from ' Barkhouse et al.l ( l2007h and comprised a set of 
57 Abell clusters over a range of richnesses. 



3 PREDICTIONS FROM THE GALFORM 
SEMI-ANALYTIC MODELS 

We draw samples from the Millennium sim ulation using the semi- 
analytic prescription of iBower et al] ( |2006|) . hereafter B06. Galax- 
ies are selected based on their halo mass (>10^'' M© for clusters) 
and their observed colours and magnitudes in this prescription. 
Red-sequence galaxies were selected by fitting a Gaussian to the 
red peak of the [B — V) colour distribution and locating the local 
minimum on the blue side to isolate the extent of the red-sequence. 
In order to capture the uncertainties in sampling the simulation in 
the same way as the observations, we vary the slope and normali- 
sation of this red-sequence limit (by 15% and 3% respectively) as 
well as the magnitude limits used to define DGR (0.1 mags each 
limit) and Monte-Carlo the results. The points for clusters and field 
in three redshift bins (0.0, 0.5, 1.0) are shown as asterisks and tri- 
angles respectively in Fig. 1. Statistical errors are negligible, so the 
error bars reflect this systematic variation. 

When confronted with the observational data, the semi- 
analytic model fails to reproduce either the redshift evolution of 
the DGR or the relative difference between clusters and the field 



at a given redshift. The model predicts approximately the correct 
DGR for clusters in the local universe, but this DGR value then re- 
mains constant out to z~ 1 . These results are insensitive to the exact 
halo mass chosen to define the clusters, as can be inferred from the 
similarity between the cluster and field points at each redshift. i.e., 
the DGR is not a strona function of halo mass in the model. 



4 DISCUSSION & CONCLUSIONS 

Previous works have suggested that there appears to be little or no 
evolution i n the bright end of the RSLF, both in c lusters and the 
field (e.g., iGilbank et alj|2007l : IScarlata et al.ll2007h . If this is the 
case, the evolution of the DGR can be attributed to the build-up 
of faint galaxies on the red-sequence with increasing cosmic time. 
This implies that star-formation is terminated in giant galaxies first 
and later in dwarf galaxies, in agreement with the "down-sizing" 
picture of star-formation. Furthermore, the lower value of DGR in 
the field compared with clusters implies that star-formation was 
completed in cluster galaxies before field galaxies. 

Performing a similar analysis on the B06 semi-analytic model 
predicts a non-evolving DGR out to z=l, with little difference be- 
tween clusters and the field. Such a result suggests that faint galax- 
ies move onto the red-sequence too efficiently in the model, with 
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this build-up of the faint end of the red-sequence already having 
been completed by the earliest epoch considered here (z~l). 

The methods available in the models for the quenching of star- 
formation include: feedback from supemovae or AGN and hot halo 
gas removal. This latter process only applies to galaxies which 
are satellites (i.e. not the most massive, central object) in their 
dark matter halo. The relative importance of the various quench- 
ing mechanisms are summarised in fig. 4 of B06. At the bright 
end of the red-sequence, most of the galaxies are the central ob- 
ject in their dark matter halo and thus feedback from AGN or su- 
pernova winds are required to turn off star-formation. At the faint 
end (in the dwarf regime), red-sequence galaxies are predominantly 
satellites. We have verified that this is still the case for all the red- 
shifts we_conside£hereJn most galaxy formation models (e.g., 
B06, iDe Lucia et d]|200Q) . when any galaxy enters a more mas- 
sive dark matter halo and becomes a satellite, it instantly has all of 
its ho t gas remo ved in a process generally referred to as 'strangula- 
tion' jLarson et al. 1980; Balogh et al. 2000). Thus, once it has con- 
sumed the cold gas already present in its disk, its fuel can no longer 
be replenished by hot gas from the halo and its star-formation is ter- 
mi nated. This prescription may be too efficient. Recent simulations 
by [McCarthy et all j2007h found that the removal of the initial hot 
galactic halo gas does not occur instantaneously, but likely takes ~ 
1 Gyr to remove the bulk of it, and that as mu ch as 30% can still 
remain after 10 Gyr. In addition, the models of lOekel & BimboimI 
( l2006h predict the existence of a critical halo mass (- 10^^ M0)for 
haloes, below which strangulation may not occur. Thus, we expect 
additional parameters involving strangulation efficiency/timescale 
and/or a threshold halo mass for strangulation to be important in 
tuning the models to reproduce the observations presented in Fig. 1 . 

To illustrate the effect of a threshold halo mass for strangu- 
lation, we use the following toy model to predict the evolution 
of the field DGR. Field surveys do not just probe isolated (i.e. 
non-group) galaxies. Indeed, they are dominated by groups, espe- 
cially toward lower redshift as cosmic structure grows. Thus it is 
group-scale physics which is important in determining the proper- 
ties of galaxies in field surveys. We assume further that galaxies in 
haloes below some (group-size) threshold mass, Alth, have no red- 
sequence satellites, i.e. DGR — 0. We assume that haloes above 
Mth have a DGR given by that of the observed cluster population 
at the same redshift (from the best-fit cluster line in Fig. 1). We 
consider Mth = (10^^, 10^^'^, 10^^) Mq and measure the frac- 
tional contribution of haloes in the appropriate mass range from 
the Millennium Simulation. We overplot the expected evolution of 
the field DGR for the different threshold masses in Fig. 1. It can be 
seen that both the Mth = 10^^ and 10^^'^ Mq models give good 
agreement with the observed field DGR. In order to reproduce the 
observations so well, our model requires both an increasing con- 
tribution to the field by groups due to growth of structure and an 
intrinsic evolution of the DGR in haloes above Mth- If either of 
these factors is removed (e.g. if a constant group fraction at all red- 
shifts or a constant cluster DGR at all redshifts is assumed) then 
the predicted evolution of the field DGR is too flat with respect to 
the observations. If groups do not exhibit exactly the same DGR 
as clusters, but instead are intermediate between clusters and the 
field (this is currently an open question), then this would effectively 
shift our predicted lines down and a lower Alth would then repro- 
duce the observations. Of course, a scenario in which the timescale 
for the removal of hot gas is much longer may fit the observations 
equally well and this or some alternative combination of threshold 
mass/timescale cannot be ruled out. Indeed, the observed evolution 



of the cluster DGR we have had to include in our model may reflect 
the effect of a strangulation timescale. 

Disentangling the effects of a timescale and threshold mass 
for the quenching of satellite galaxies is a complex observational 
problem. The colour difference between central and satellite red- 
sequence galaxies may be a useful test of models. For example, 
in the B06 prescription, it is predicted that satellite red-sequence 
galaxies are significantly redder than central red-sequence galaxies. 
,van den Bosch et al. (2007) recently measured the average colours 
of satellite and central galaxies in the local universe from an ap- 
proximately halo mass-selected sample and found that the former 
were indeed on average redder. They also showed that measuring 
halo masses from statistical samples for '--^lO^^ Mq systems is not 
possible with high completeness and thus, if there is a threshold 
mass for strangulation around this limit, suc h studies will requir e 
detailed targeted observations of groups (e.g. JWilman et all2005h . 
An alternative approach to get at the timescale for strangulation is 
the study of the star formatio n histories of indiv idual galaxies un- 
dergoing transformation (e.g.. lMoran et al.l2007l) . It is possible that 
the observations in Fig. 1, coupled with new strangulation prescrip- 
tions in the semi-analytic models may improve constraints on the 
likely values of the threshold mass and timescale for strangulation. 
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